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ABSTRACT 


An anisotropic elasticity crack tip stress analysis Is Implemented 
using three crack extension direction criteria (the normal stress ratio, 
the tensor polynomial and the strain energy density) to predict the 
direction of crack extension In unidirectional off-axis graphite- 
epoxy. The theoretical predictions of crack extension direction are 
then compared with experimental results for 15° off-axis tensile coupons 
with center cracks. Specimens of various aspect ratios and crack orien- 
tations are analyzed. It Is shown that only the normal stress ratio 
criterion predicts the correct direction of crack growth. 
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INTRODUCTION 


A fundamental problem In predicting the failure of laminated com- 
posite materials Is prediction of the direction of crack growth In the 
Individual laminae, and the laminate. The Influence of the direction of 
crack growth on the failure response of the laminate is shown In Fig. 1 
(1). The clustered (^"^s 9 ra P h1te - e P° x y laminates fall In a pure 
matrix mode (delamlnatlon and either Intralaminar matrix cracking or 
fiber matrix debonding). In contrast, the alternating [(+e/-e) 2 ) s lami- 
nates exhibit fiber breakage In half of the plies, and either matrix 
cracking or fiber matrix debonding In the others; there Is no delamlna- 
tlon. The mode of failure has a significant effect on the strength of 
the laminate. The strength of 10° and 30° alternating laminates Is, for 
example, 30 and 50 percent greater, respectively, than the strength of 
clustered 10" and 30" laminates (11. Hence, understanding the parame- 
ters that affect laminate failure, particularly those Influencing the 
direction of crack growth In the lamina and between laminae, Is of 
critical Importance In predicting the fracture response of laminates. 

Predicting the direction of crack extension In laminates Is a very 
complex three-dimensional problem. Since the lamina Is the basic build- 
ing block of the laminate. Its behavior must be fully understood as a 
stepping stone toward understanding the behavior of the laminate. This 
study was undertaken to assess more critically the applicability of 
three criteria which have been presented In the literature for predict- 
ing the direction of crack growth In unidirectional fibrous composites 
[2-4], 
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CRACK EXTENSION DIRECTION CRITERIA 
Three phenomenological criteria for predicting the direction of 
crack extension In homogeneous, anisotropic materials are the normal 
stress ratio criterion (21, the tensor polynomial criterion [3| and the 
strain energy density criterion [ 4] . These criteria can be used to 
predict the load at failure and the direction of crack extension. The 
crack tip coordinate system used In this analysis Is shown In Fig. 2. 

Normal Stress Ratio Criterion 

Buczek and Herakovlch [2| have hypothesized the normal stress ratio 
criterion as a crack growth direction criterion. The model assumes that 
the direction of crack extension corresponds to the direction of the 
maximum value of the normal stress ratio R(r o ,$) where 

R(r„,*) - (1) 

0 T 

<H 

In the expression for R(r a,, corresponds to the normal stress 

acting on the radial plane defined by $, and at a given distance, r Q , 
from the crack tip. T. . is the tensile strength on the * plane. 

Since the tensile strength on an arbitrary plane Is difficult, if 
not impossible, to measure, T Is defined in a manner consistent with 

99 

the tests that can be performed. To meet this requirement, a mathe- 
matical definition of T.. must satisfy the following cor.dltlons: 

99 

(1) for an Isotropic material, T xx must be Independent of <t> . 

<9V 

(2) for crack growth parallel to the fibers, T^ must equal the trans- 
verse tensile strength Yy. 
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(3) for crack growth perpendicular to the fibers, T. , must equal the 
longitudinal tensile strength Xy. 

A definition satisfying these conditions is; 

* Xy sin 2 a + Yy cos 2 b (2) 

where s is the angle from the plane of Interest to the fiber direction. 

Tensor Polynomial Criterion 

Tsai and Wu (5] first presented the tensor polynomial criterion as 
an anisotropic failure criterion. This criterion is based on the 
existence of a failure surface In stress space of the form: 

f (3) 

where Fj and F^j are strength tensors of second and fourth order, 
and oj Is the contracted form of the stress tensor. Expressions for 
and are given In Table 1. 

In application of the tensor polynomial to fracture problems (3), 
the assumed direction of crack extension Is the radial direction of 
maximum f(o^). The stress components are those determined by a 

continuum mechanics-based stress analysis, and must be evaluated at a 
finite distance, r Q , from the crack tip. 

Strain Energy Density Criterion 

The strain energy density criterion Is based on variations in the 
energy stored along the periphery of a core region surrounding the 
crack. S 1 h presents the criterion for Isotropic fracture In [6] and a 
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modified form for application to anisotropic fracture In [4|. 
The strain energy density factor, S, Is defined as: 


3W „ S 

17 7 


( 4 ) 


where |rf Is the strain energy density function and r Is the distance 
from the crack tip. Since the strain energy density function can be 
expressed in terms of the crack tip stresses and strains for plane 
stress as: 


W ’ 5 ( Vx + Vy + 'xyV <5) 

an expression for the strain energy density factor, S, can be obtained 
by substitution. The resulting expression Is: 

5 ’ 1 <Vx + Vy + Vxy’ (6) 

The fundamental hypothesis of Sih (4j for unstable crack growth Is that 
crack Initiation takes place In the radial direction corresponding to a 
minimum value of the strain energy density factor, i.e., 

2 

-I = 0 and -2-f >0 at <t> = (7) 

3 * a/ c 

Sih cautions that for small values of r, a continuum mechanics- 
based crack tip stress analysis is Invalid. Hence, the strain energy 
factor should be evaluated at a finite distance, r QI from the crack tip, 
where r Q Is of the same order of magnitude as the crack tip curvature. 
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ANISOTROPIC ELASTICITY ANALYSIS OF CRACK TIP STRESS FIELDS 
The stress analysis of an Infinite homogeneous anisotropic plate 
containing a center crack can be directly related to a homogeneous 
anisotropic p’ate with an elliptic hole. By reducing the minor axis 
dimension to zero and evaluating the stress potential functions In the 
neighborhood of the crack tip, Lekhnl t sk 1 1 1 s complex variable solution 
(7) for an elliptic hole In an anisotropic plate can be adapted to solve 
anisotropic fracture problems. Wu presents a detailed descn tlon of 
this procedure In [8J, along with equations describing the crack tip 
stresses for an Infinite homogeneous anisotropic center-cracked plate. 
The problem under consideration Is shown In Fig. 3. 

The governing partial differential equation for this problem In 
terms of the Airy's stress function U Is: 


a 4 U 2A 26 a 4 ll (2A 12 + A 66 ) a 4 U 2A 16 3 4 U Ajj 3 4 U 

- _ 3 • + p o “ J + 4 

ax' A 22 ax J ay A 22 ax^ay^ A ?2 axay J A 22 ay 


( 8 ) 


where are components of the compliance tensor for plane stress or 
plane strain, depending on the analysis desired. 

Assuming U = e x+sy , the characteristic equation for (8) takes the 

form: 


AnS 4 - 2A lg S 3 + (2A 12 + A 66 )S 2 - 2A 26 S + A ?2 = 0 (9) 

The roots of the characteristic equation, and S 2 , (and their conju- 
gates) are complex, and are functions of the material properties and the 
orientation of the crack relative to the principal material direction. 
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Assuming Sj + S 2 , evaluation of the complex potential functions 
near the crack tip yields expressions for the stress and displacement 
distributions of the form: 


o" /a S,S P S ? 

Re ( — L?- |-i 

/5f (Sj - s z ) *![ 


t /a 

+ Re{- 


1 S 2 s 2 

— [---II 


/Tr (Sj - Sj) q 


♦i 


o- /a 1 ,S. S« t /a 1 » i 

« y Re ( [-5- - ^|( * Re [ [- 1 - - - 1 -]) 


✓?r 


(Sj - S 2 ) n), 


✓2? 


(Sj - S 2 ) * 2 


ot> “ 

o /a 


S,S„ 1 1 

Re { — L t — | -]) * 


00 /** 
t /a 


1 S, S, 

Re ( [-*--£]) 


/2r - S 2 ) ^ ^ {S 1 " S 2^ *1 *2 


u » o“ /2ar Re { -^g - S~)^ S 1 P 2^2 " S 2 p l' i 'l^ 

+ t“ /I5r Re ( ( $i \ $ 2 y lP 2 »jj ' P^i 11 


( 10 ) 


/2ar Re {(s i - S 2 )f S l q 2 ll/ 2 ' S 2 q l il l^ 

+ t” /2ar Re ~ s g ) t q 2^2 “ q l*l ^ 


= cos* + S^sln* 


P 1 3 A 11 S 1 + A 12“ A 16 S 1 


A 12 S 1 + A 22“ A 26 S 1 


t|> 2 = cos* + S 2 s1n* 

p 2 = a u s 2 2 + a 12 - a 16 s 2 

_ A 12 S 2 2+ A 22" A 26 S 2 
q 2 = ; 


where 
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As In Ihe Isotropic case, the crack tip stresses exhibit a 
singularity of l//r. However, the magnitude of the stresses Is not 
simply a function of the stress Intensity factors. The quantities Sj 
and S 2 also affect the magnitude of the stresses. This Is an Important 
difference between anisotropic and Isotropic fracture. In anisotropic 
fracture, the magnitude of the crack tip stresses Is a function of not 
only the applied load, specimen geometry and crack length, but also the 
material properties and the orientation of the crack relative to the 
principal material direction. Application of this solution to the 
analysis of unidirectional composites with crack orientations other than 
parallel to the X-axis Is presented In Gregory and Herakovlch (9). 

EXPERIMENTAL PROGRAM 

In order to test the ability of the theoretical models to predict 
the direction of crack extension, tests were performed on center-cracked 
specimens of 16 ply unidirectional AS4/3501-6 graphite/epoxy. Material 
properties for AS4/3501-6 graphite/epoxy are given in Table 2. The 
experimental Investigation consisted of a series of 15° off-axis tensile 
tests with rigid end constraints. Coupons of various aspect ratios with 
pre-machlned cracks oriented perpendicular to the loading direction or 
perpendicular to the fibers were tested. The specimen configurations 
are Illustrated In Fig. 4; two coupons were tested for each 
configuration. 

The test specimens were 1.00 Inch (25.4 mm) wide and contained 
center cracks of 0.20 Inches (5.08 mm). The aspect ratios (gauge length 
to width) of the specimens were 8, 4, and 1. Since off-axis specimens 
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with rigid end constraints experience Increased shear loading as the 
aspect ratio Is reduced (10), specimens of various aspect ratios were 
tested to vary the biaxial stress state In the region of the cracfc. 

Strain gauge rosettes were attached to each specimen at a point far 
from the crack. During each test, the direction of crack extension, the 
load at crack Initiation* and the load at failure were observed and 
measured. To facilitate visual observation of crack growth, load was 
Incrementally applied at a crosshead speed of 40 microstrain per second. 

CORRELATION OF THEORETICAL MODELS WITH EXPERIMENTAL RESULTS 

In all the experiments performed, crack extension occurred parallel 
to the fibers, with no observable fiber breakage. A broken specimen for 
each combination of crack orientation and aspect ratio Is shown In Fig. 
5. 

The anisotropic elasticity solution was used to model the experi- 
mental procedures. To better approximate the far field stresses acting 
on the crack, the Pagano and Halpln solution (10) for the state of 
stress In an off-axis tensile coupon with rigid end constraints was 
Incorporated. Though the Pagano and Halpln solution does not account 
for the presence of a crack. It does account for the high stress 
gradients and shear stress generated by the fixed ends. The far field 
stresses used as Input for the anisotropic elasticity solution are the 
stresses generated by the Pagano and Halpln solution at a point 
corresponding to the crack tip location. The far field stresses, 
predicted direction of crack growth and the experimentally observed 
direction of crack extension are given In Table 3. 
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From Table 3, It Is apparent that only the normal stress ratio 
criterion predicts the correct direction of crack extension. The other 
crack extension direction criteria show no correlation with the experi- 
mental results. Distributions In the crack extension criteria a? a 
function of 4 >, for test case A (crack perpendicular to the loading 
axis), with aspect ratio 1 are shown In Figs. 6-8. Distributions In the 
normal stress ratio as a function of for aspect ratios 1, 4, and 8 of 
test case B, are shown In Fig. 9-11. 

Analysis of Figs. 6, 9, 10 and 11 and Table 3 yields an interesting 
characteristic of the normal stress ratio criterion. The theoretical 
predictions of crack extension direction differ slightly from the exper- 
imentally observed values. There is, however, a strong peak In the 
distribution of the normal stress ratio as a function of $ In the actual 
direction of crack growth. This fact Is very Important. The normal 
stress ratio may not have the accuracy to predict correctly the 
direction of crack extension to within one degree. When observed graph- 
ically, however, the normal stress ratio represents the direction of 
crack extension exceptionally well. 

The normal stress ratio criterion predicted the correct direction 
of crack extension for every test analyzed, except for test case B 
(crack perpendicular to the fibers) with aspect ratio 1. For this 
problem, the normal stress ratio correctly predicts crack extension 
parallel to the fibers, however the predicted direction of extension is 
180° out of phase with the observed direction. Analysis of Figs. 9-11 
reveals that for test case B, there are two peaks In the normal stress 
ratio. The first peak, near $ = -90°, predicts crack extension parallel 
to the fibers, toward the center of the coupon. The second peak, 
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near $ » +90°, Implies crack extension parallel to the fibers, toward 
the free edge. For aspect ratios 4 and 8, the second peak Is the 
maximum value. This Is not true for aspect ratio 1, which has a maximum 
value at the first peak. 

The discrepancy In the normal stress ratio criterion for test case 
B with aspect ratio 1 does not necessarily compromise the validity of 
the criterion. The stress gradients in a rigidly constrained tensile 
test of aspect ratio 1 are very high. One can also question the ade- 
quacy of the crack tip stress analysis; 1.e.» using Pagano and Hal pin 
tensile coupon stresses as far field stresses In the anisotropic elas- 
ticity solution. The discrepancy Is noted and further research Is 
required. 


CONCLUSIONS 

This study was concerned with the development of a model to predict 
the direction of crack extension In unidirectional composite 
materials. An anisotropic elasticity solution, In conjunction with 
Pagano and Halpln's solution for stresses In a fixed end tensile test, 
was used to calculate the crack tip stress field In a center-cracked 
off-axis tensile coupon. Three crack extension direction criteria, the 
normal stress ratio criterion, the tensor polynomial criterion and the 
strain energy density criterion, were then incorporated into the model 
to predict the direction of crack extension. 

Comparison of the predicted direction of crack extension with 
experimentally observed crack growth. Indicates that only the normal 
stress ratio criterion consistently predicts the correct direction of 
crack extension. 
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Table 1 

Relationships for Strength Tensors In Terms of Measured Strengths 


Fi - (1/X T + 1/X C ) 
F 2 » (1 /Y t + 1/Y C ) 

Fg = 0.0 

F U = -1/(X T X C ) 

F 22 = -i/(v t y c ) 
f 6 6 ■ 1/(S 2 ) 


Table 2 

Lamina Properties of AS4/3501-6 Graphite/Epoxy 

E t * 21.6 MSI (148.9 GPa) 

E 2 = 1.96 MSI (13.5 GPa) 

G 12 = 0.83 MSI (5.7 GPa) 

X T = 282 KSI (1.94 GPa) 

X c = -282 KSI (-1.94 GPa) 

Y t = 10 KSI (68.9 MPa) 

Y c = -10 KSI (-68.9 MPa) 

S = 14.2 KSI (97.9 MPa) 
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(v) The far field loads correspond to the crack tip loading at 136 strain. 
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FIG. 5 FAILED CENTER CRACKED SPECIMENS 
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FIG. 7 TENSOR POLYNOMIAL NEAR A CRACK TIP, AR-|, CASE A 
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FIG, 9 NORMAL STRESS RATIO NEAR A CRACK TIP, AR=1, CASE B 
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FIG. 10 NORMAL STRESS RATIO NEAR A CRACK TIP, AR=4, CASE B 
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FIG. 11 NORMAL STRESS RATIO NEAR A CRACK TIP, AR=8, CASE B 
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